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ABSTRACT 

We present the results of Magellan/MMIRS and Keck/NIRSPEC spectroscopy for five Lya emitters 
(LAEs) at z ~ 2.2 for which high-resolution FUV spectra from Magellan/MagE are available. We 
detect nebular emission lines including Ha on the individual basis and low-ionization interstellar (LIS) 
absorption lines in a stacked FUV spectrum, and measure average offset velocities of the Lya line, 
Ai>Lya, and LIS absorption lines, Ati a b s , with respect to the systemic velocity defined by the nebular 
lines. For a sample of eight z ~ 2 — 3 LAEs without AGN from our study and the literature, we obtain 
A^Lya = L75 ± 35 km s _1 , which is significantly smaller than that of Lyman-break Galaxies (LBGs), 
A^Lya — 400 km s _1 . The stacked FUV spectrum gives Au a b s = —179 ± 73 km s _1 , comparable to 
that of LBGs. These positive AuLya and negative Au a b s suggest that LAEs also have outflows. In 
contrast to LBGs, however, the LAEs' Av-^ ya is as small as |Ai> a b s |, suggesting low neutral hydrogen 
column densities. Such a low column density with a small number of resonant scattering may cause 
the observed strong Lya emission of LAEs. We find an anti-correlation between Lya equivalent width 
(EW) and Av^ ya in a compilation of LAE and LBG samples. Although its physical origin is not clear, 
this anti-correlation result appears to challenge the hypothesis that a strong outflow, by means of a 
reduced number of resonant scattering, produces a large EW. If LAEs at z > 6 have similarly small 
Ai>Lya values, constraints on the reionization history derived from the Lya transmissivity may need 
to be revised. 

Subject headings: cosmology: observations — galaxies: formation — galaxies: evolution — galaxies: 
high-redshift — 



1. INTRODUCTION 

Lya emitters (LAEs) are objects with a large rest- 
frame Lya equivalent width, EW(Lya) > 20-30 A. This 
population is usually selected using a narrow band filter 
for Lya emission combined with a broad band filter which 
measures continuum emission around Lya. Very re- 
cently, nearby LAEs down t o z ~ 0.2 have also been stud- 
ied with GALEX data (e.g. JCowie et al.ll201lD . Previous 
studies have revealed that most (high-z) LAEs are young, 
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low-mass galaxies with smal l dust extinction, while some 
are old, massive, and dusty (|Ono et a l. 2010a). Morpho- 
logical studies have revealed that LAE s are smaller than 
typic al high-z star forming galaxies (|Bond et al.l 120091 
12010th Thus, LAEs are among the building block candi- 
dates in the A Cold Dark Matter (CDM) model, where 
smaller and le ss massive galaxies merge to be larger and 
massive ones (jRauch et al.l[2008l ). 

Gas exchanges between galaxies and the ambient in- 
tergalactic medium (IGM), i.e., outflows and inflows, are 
thought to play important roles in galaxy evolution. Out- 
flows are driven by supernovae (S Ne), stellar winds from 
massive stars , and AGN activity dHeckman et all 119901 : 
IMurrav et all 120051: IChoi fc Nagamind 120 111) . Galaxies 
which lost cold gas via outflows may experience a re- 
duction or even termination of subsequent star forma- 
tion. In contrast, cold gas supplied by inflows may in- 
crease the star formation, especially when they comes in 
the form o f dense, filamentar y gas streams ('cold accre- 
tion'; e.g. IDekel et al.l 120091 ). Gas exchanges will also 
affect the chemical evolution of galaxies and the IGM, 
and thus have played important roles in establis hing the 
mass - metalliciy relation ac r oss cosmic time ( e.g., iLarsonl 
[l^lTmnonti et al.|[2CjollErb et al.ll2006al) . 

Outflows have been found in nearby sta rburst galaxie s 
(jHeckman et al.l 119901 ) . nearby ULIRGs (iMartinl I2005D 
Lyman-break galaxies (LB Gs) at z ~ 2 — 3 (jPettini et al.l 
12001 iShaplev et al l[200l . and BX/BM galaxies at z ~ 2 
(|Steidel et al.ll2010l ). These studies have made use of the 
fact that FUV low-ionization interstellar (LIS) absorp- 
tion lines, which are generated when continuum photons 
encounter the outflowing gas, are blue-shifted with re- 
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spect to the systemic redshift measured by nebular emis- 
sion lines such as Ha originated from Hn regions in the 
galaxy. ^ 

Examining the incidence of outflows in LAEs is of inter- 
est, because LAEs are generally less massive (< IQ 9 M^; 
Nilsson et al.1 120071: iFmkelstein et all 120081 : lOno et alJ 
2010allbl : INakaiima et al.ll2012al) than, e.g., LB Gs (10 10 - 
10 11 M^: [ShaDlev et al.l[200l lErb et al"1l2006bl) and BzK- 
selected galaxies (10 10 -10 n M^; e.g.. lDaddi e"t~al][20p: 
lYuma et al.ll2012f ) , and thus have shallower gravitational 
potentials. Some authors have argued the importance of 
outflows fro m less massiv e galaxies in chemical enrich- 
ment (e.g.. iLarsonl I1974T) . However, FUV continua of 
LAEs are too faint for LIS absorption lines to be reliably 
measured with current facilities. 

Lya is also used to probe the gas kinematics. The 
Lya line is known to have complex profiles caused by 
its resonant nature. Many theoretical and observational 
studies have shown that outflowing gas leads to a red- 
shift ed Lya line with respect to the systemi c redshift 
(e.g.. lVerhamme et al.|[2006HSteidel et al.ll2010[) , giving a 
positive Lya offset velocity, AvLya ■ In the case of an out- 
flow, back-scattered (i.e., redshifted) Lya photons have 
more chance of escape because they drop out of resonance 
with the foreground gas. However, it is very difficult to 
obtain high-resolution nebular emission spectra for LAEs 
because they are faint compared to strong sky emission. 

Prior to our study, only four LAEs have both high- 
resolution Lya and n ebular line sp ectra: two from 
iMcLinden et al.1 ()2011l ) and two from iFinkelstein et all 
(|2011| ). McLinden et al.'s objects have Avh ya = 125 ± 
17.3 and 342 ± 18.3 km s _1 , respectively, and Finkelstein 
et al.'s objects have A^Lya = 288 ± 37 (photometric) ± 
42 (systematic) and 189 ± 35 ± 18 km s _1 , respectively 
(see also Chonis et al., in preparation). The fact that all 
four have Av-^ ya > suggests that outflows are common 
in LAEs, but due to the small sample sizes obtained so 
far, there has not been a statistical discussion of the gas 
motions of LAEs. 

Answering the fundamental question as to why LAEs 
have strong Lya emission is crucial for understanding the 
physical nature of LAEs. Some studies have shown that 
outflows facilitate the escap e of Lya photons from galax- 
ies (e.g.. iKunth et al.lll998fi as they reduce the number 
of re sonant scattering. Indeed, iVerhamme et al.l ((20061 
l2008f ) have carried out Lya radiative transfer simula- 
tions, and claimed that expanding shell models can ac- 
count for observed Lya spectral profiles of LBGs. How- 
ever, while certain aspects of the Lya radiative transfer 
mechanism of high-z galaxies can be understood in the 
context of outflows, the reason for the strong Lya emis- 
sion of LAEs is still an open question. 

In order to address these questions, we are conduct- 
ing near infrared spectroscopy for optically confirmed 
LAEs at z ~ 2.2 promising for detecting nebular emis- 
sion lines. High-resolution spectra of nebular emission 
lines are essential for measuring the offset velocities of 
the Lya line, Ai>Ly a , and of LIS absorption lines, Au a b s , 
with respect to the systemic velocity. The redshift of 
z ~ 2.2 is favored since we can simultaneously observe 
Lya, LIS absorption lines (e.g., [Sin]), [Oi], and [Cn]), 
and optical nebular lines (e.g., [On], H/3, [Oiii], and Ha) 
from the ground. Furthermore, we can compare the kine- 



matics and Lya radiative transfer of LAEs with those of 
brighter and more massive galaxies at similar redshifts, 
e.g., LBGs, obtained by previous studies. 

We successfully detect for several LAEs nebular emis- 
sion lines on the individual basis and LIS absorption 
lines in a stacked FUV spectrum, to measure At>L yQ 
and Aw a bs- We also derive physical quantities such as 
EW(Lya), Lya escape fraction (f^/ c a ), star formation 
rate (SFR), stellar mass, and the size of the star- forming 
region from the spectral data and photometric data (SED 
fitting). Using these quantities, we discuss the gas kine- 
matics of LAEs and give implications on the physical 
origin of the strong Lya emission. 

This paper is organized as follows. Our near infrared 
spectroscopy is described in $U After performing SED 
fitting in Sj3j we derive At>Ly Q , At> a b s , and several ob- 
servational quantities for our LAEs in fj4j A discus- 
sion in the context of outflow and Lya radiative trans- 
fer is given f|5l followed by conclusions in ij6] Through- 
out this paper, mag nitudes are given in the AB system 
(|Oke fc Gunnlll983l) . and we assume a ACDM cosmology 
with fl m = 0.3, fl A = 0.7 and H = 70 km s" 1 Mpc" 1 . 

2. SPECTROSCOPIC DATA 

2.1. Targets of our NIR Spectroscopy 

Our targets for near infrared spectroscopy are selected 
from samples of z ~ 2.2 LAEs in the COSMOS and 
the Chandra Deep Field South (hereafter CDFS) con- 
structed by K. Na kajima et al. (in prepa ration) in the 
same manner as m iNakajima et al.1 (|2012af) . These LAE 
samples are based on narrow-band (NB387) imaging with 
Subaru/Suprime-Cam, supplemented by public broad- 
band data. NB387 is our custom-made filter with a 
central wa velength and FWHM o f 3870 A and 94 A, re- 
spectively (Nakaii ma et al.ll2012al ). The LAEs have been 
selected by imposing the following color criteria: 

u* - NB387 > 0.5 & B - NB387 > 0.2 (COSMOS)(l) 
U - NB387 > 0.8 & B - NB387 > 0.2 (CDFS). (2) 

The COSMOS (CDFS) sample contains 619 (1,108) 
LAEs with EW>30A down to NB387 = 26.1 (26.4). 
Among them, two objects in the CDFS sample (ODFS- 
3865, CDFS-6482) and three in the COSMOS sample 
(COSMOS-13636, COSMOS-30679, COSMOS-43982) 
have a high-quality Magellan/MagE spectrum of Lya 
(M. Rauch et al., in preparation). These five objects 
are our targets of NIR spectroscopy. They have typical 
NB387 excesses among the whole LAE sample in each 
field while having relatively bright NB387 magnitudes. 

2.2. Near Infrared Spectroscopy 

We observed the two CDFS objects on 2010 Octo- 
ber 21 with Magellan/MMIRS using the HK grism cov- 
ering 1.254 - 2.45/zm. The total exposure time was 
10800s for each object. The slit width was 0."5 result- 
ing in R = A/AA ~ 1120. A two-point dither pattern 
(A1,B1,A2,B2,A3,B3,...) was adopted. The A0V stan- 
dard star HIP-16904 was also observed. The sky was 
clear through our observation run, with seeing sizes of 
0."5 - 0."9. 

The three COSMOS objects were observed on 2011 
February 10 and 11 with Keck-II/NIRSPEC. COSMOS- 
30679 was observed with NIRSPEC-3 (J band; 1.15 - 
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1.36/im), NIRSPEC-5 (H band; 1.48 - 1.76/mi), and 
NIRSPEC-6 (X band; 2.2 - 2.43/im) filters in the low- 
resolution mode, while COSMOS-13636 and COSMOS- 
43982 were observed with the K band alone. Further- 
more, we observed CDFS-3865 with the J band target- 
ing [On] AA3726, 3729. Total exposure times are shown 
in Table HJ The slit width was 0."76 for all three ob- 
jects, corresponding to R ~ 1500 for all grisms. A two- 
point dither pattern was adopted. We simultaneously 
obs erved reference stars, which we used for b lind offsets 
(see lFinkelstein et al.ll20TTt lYang et al.ll2011| ) because of 
the faintness of our targets. The AOV standard star HIP- 
13917 was also observed. The sky was clear in our observ- 
ing nights, with seeing sizes of 0."6 - 0."9. A summary 
of the spectroscopic observations is given in Table [T] 

2.3. Data Reduction 

We reduced the MMIRS data using IRAF tasks and 
the COSMOS package which is the standard reduction 
pipeline for Magellan/IMACS. The MMIRS detector is 
read out non-destructively during a single exposure, and 
individual read-outs are stored as separate extensions in 
a FITS file through which we know whether and when a 
particular pixel is saturated. Bias subtraction and flat 
fielding were processed for each read-out using IRAF 
mscred package to treat data of this format. Then we 
ran mmf ixen package which takes advantage of its sam- 
pling. This essentially fits a line to different values for 
a given pixel in each readout, and outputs the slope of 
this linear fit to the final collapsed image. Wavelength 
calibration and distortion correction were processed for 
each frame using the COSMOS package. Although we 
obtained arc lamp calibration images, we used OH lines 
for wavelength calibration. We then performed the fol- 
lowing operation to remove sky background: CI = Bl - 
(Al + A2)/2, C2 = B2 - (A2 + A3)/2. After this op- 
eration, we ran subsky in the COSMOS package E3 on 
each frame to remove residual sky lines. Resultant frames 
(CI, C2,...) are then stacked to have a final 2D frame us- 
ing sumspec-2d and extract-2dspec in the COSMOS 
package. ID spectrum extraction was carried out using 
apall in IRAF. The telluric absorption correction and 
flux calibration were conducted using the standard star 
frames. The flux-calibrated ID spectra of the two objects 
are shown in Figure [TJ 

We reduced the NIRSPEC data using mainly IRAF 
tasks. Details o f the red uction procedure is described in 
iNakaiima et al.l (j2012b[ K The flux-calibrated ID spectra 
of the three objects are shown in Figure [5J 

2.4. Emission Line Detections and Measurements 

We determine a line to be detected, if there exists an 
emission line above the 3<r sky noise around the wave- 
length expected from the Lya redshift, where sky noise 
is calculated from the spectrum within 100 A from the 
line wavelength. Table [2] summarizes the line detections. 
For each detected line, we fit a Gaussian function to cal- 
culate its central wavelength, FWHM, and flux. The cen- 
tral wavelength is then converted into the vacuum wave- 
length after correction for the heliocentric motion of the 
earth EE and is adopted for the systemic redshift of the 
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Fig. 1.— Reduced 1D+2D spectra of CDFS-3865 (top) and 
CDFS-6482 (bottom). 
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11 This uses the Kelson procedure (cf, Kelson 2003). 

12 http: / /fuse. pha.jhu.edu/support/tools/vlsr. html 
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Fig. 2.— Reduced 1D+2D spectra of COSMOS-13636 (top), 
COSMOS-30679 (middle), and COSMOS-43982 (bottom). H-band 
data were taken only for COSMOS-30679. 

object. For the four objects with multiple line detections, 
we confirm that the redshifts of the lines agree within la 
errors, and adopt their weighted mean for the systemic 
redshift. [On] is not used in the redshift measurement 
because the 3726 A/3729 A doublet is not resolved in our 
spectra. The H/3 line for CDFS-3865 is not used, either, 
because of the low S/N ratio. We also derive la upper 
limits of fluxes for all undetected lines. The emission line 
measurements are summarized in Table [2j 
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TABLE 1 

Summary of our observations 



Object 


a(J2000) 


<5(J2000) 


z(Lya) 


L(Lya) 


Dates 


^cxp 










(10 42 erg s- 1 ) 




00 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 



CDFS-3865 03:32:32.31 -28:00:52.20 2.17507I ( ~ 29.8 ± 4.9 2010 Oct 21 5100 (J), 10800 (HK) 

CDFS-6482 03:32:49.34 -27:59:52.35 2.20610t°°g° U 2 15.4 ±8.1 2010 Oct 21 10800 (HK) 

COSMOS-13636 09:59:59.38 +02:08:38.36 2.16229 ± 0.00008 11.3 ±0.5 2011 Feb 10-11 5400 (K) 

COSMOS-30679 10:00:29.81 ±02:18:49.00 2.20046 ± 0.00008 8.5 ± 0.7 2011 Feb 10-11 5400 (J), 7200 (H), 6300 (K) 

COSMOS-43982 09:59:54.39 ±02:26:29.96 2.19396 ± 0.00008 11.0 ±0.5 2011 Feb 10-11 3600 (K) 



Notes. (1) Object ID; (2), (3) Right ascension and declination; (4) Rcdshift of Lya emission; (5) Lya luminosity derived from narrow- 
and broad-band photometry; (6) Dates of observations; (7) Exposure times for the filters shown in parentheses. 



TABLE 2 
Emission Line Properties 



Object 


Line 


Arcst 


^obs 


Acorr . 


z 


S/N 


Insuturmcnt 






(A) 


(A) 


(A) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




[On] 


3727.00 


a 


a 


(7, 


6.2 


NIRSPEC 




H/3 


4861.33 


15418.9 ±4.0 


15418.3 


2.17162 ±0.00083 


3.2 


MMIRS 




[Oiii] 


4958.91 


15735.8 ±3.5 


15735.1 


2.17310 ±0.00071 


5.4 


MMIRS 


CDFS-3865 


M 


5006.84 


15887.1 ± 1.5 


15886.4 


2.17294 ±0.00030 


15.3 


MMIRS 




Ha 


6562.85 


20825.8 ± 1.7 


20824.9 


2.17315 ±0.00026 


16.0 


MMIRS 




[Nil] 


6583.45 








< 1 


MMIRS 



H,3 4861.33 - - - < 1 MMIRS 
[Oiii] 4958.91 - - - < 1 MMIRS 
CDFS-6482 [Oiii] 5006.84 16048.3 ± 3.8 16047.63 2.205144 ± 0.00076 8.8 MMIRS 
Ha 6562.85 21036.9 ± 1.9 21036.02 2.205317 ± 0.00029 5.5 MMIRS 
[Nil] 6583.45 - - - < 1 MMIRS 



COSMOS-13636 Ha 6562.85 20753.5 ± 1.0 20752.36 2.16210 ± 0.00015 7.1 NIRSPEC 
[No] 6583.45 - - - < 1 NIRSPEC 



[On] 3727.00 - a - a - a 4.0 NIRSPEC 

U/3 4861.33 - - - - b NIRSPEC 

[Om] 4958.91 - - - - b NIRSPEC 

COSMOS-30679 [Om] 5006.84 16013.8 ± 1.2 16013.62 2.19835 ± 0.00024 10.0 NIRSPEC 

Ha 6562.85 20993.7 ± 1.3 20993.48 2.19884 ± 0.00019 11.5 NIRSPEC 

[Nil] 6583.45 - - - <1 NIRSPEC 



COSMOS-43982 Ha 6562.85 20958.50 ± 1.2 20958.26 2.19347 ± 0.00018 11.4 NIRSPEC 
[Nil] 6583.45 21026.20 ±2.1 21025.96 2.19376 ± 0.00032 7.1 NIRSPEC 

Note. — The weighted mean redshifts for objects with multiple line detection are 2.17306 ± 0.00019 
(CDFS-3865), 2.20530 ± 0.00027 (CDFS-6482), 2.19865 ± 0.00015 (COSMOS-30679), and 2.19354 ± 0.00016 
(COSMOS-43982). The symbol "-" indicates no detection. (1) Object ID; (2), (3) Line name and its rest- 
frame wavelength; (4) Observed wavelength of the line; (5) Wavelength of the line corrected for the LSR 
motion; (6) Redshift; (7) Signal to noise ratio of the line detection; (8) Instrument. 

a [On] redshift is not shown because it is not reliably measured. 

b S/N upper limit is not shown because the line is contaminated by a strong OH line. 



2.5. Checking the presence of AGNs 

We examine whether our objects host an AGN in three 
ways. First, we compare the sky coordinates of the ob- 
jects with those in deep archival X-ray and radio cat- 
alogues. For CDF S , we refer to the X-ray catalogue 
given by iLuo et al.l (|2008f ). whose sensitivity limits arc 
1.9 x 10" 17 and 1.3 x 10" 16 ergs cm" 2 s" 1 for the 0.5-2.0 
and 2 - 8 keV bands, respec tively. B For CO SMOS, 
we use the X-ray catalogue bv lElvis et al.l (|2009D . whose 
sensitivity limits are 1.9 x 10~ 16 (0.5 — 2.0 keV band), 
7.3 x 10~ 16 (2 - 10 keV band), and 5.7 x 10~ 16 ergs 
cm~ 2 s^ 1 (0.5 — 10 keV b and) . We also ref e r to t he radio 
catalogue constructed by ISchinnerer et al.l (|2010D , whose 



sensitivity limits are 10 — 40/^ Jy/beam. No counter- 
part for the LAEs is found in any of the catalogues. 

Second, we look for three high ionization state lines 
typical of AGNs, CIVA1549, HeIlA1640, and CIIL1A1909, 
in the spectra, and detect none of them. 

Finally and m ost importantly, we a pply the BPT diag- 
nostic diagram ([Baldwin et al.l [l98lh to our objects, as 
shown in Figure [3] The solid curve in Figure [3] shows 
the boundar y between sta r -formi ng galaxies and AGNs 
proposed bv iKewlev et al.l (|2001| ) using photoionization 
models, while the dotted curve is the boundary empir i- 
cally defined by SDSS objects (jKauffmann et al.l 120031) . 
In both cases, star-forming galaxies fall below the curve. 



CDFS |http://www2.astro.psu.edu/users/niel/cdfs/luo2OO8-chandra-cE felC^g>M0ff IMtpd/s'4Iaa.tpMo2attdch .edu/data/COSMOS / 1 
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None of our objects has detections of all four lines nec- 
essary for the BPT diagram. Thus, as shown in Figure 
[3l any object but one logically has a possibility of being 
in the AGN regime. However, we conclude from the fol- 
lowing discussion that all but COSMOS-43982 are star- 
forming galaxies. CDFS-3865 falls below the curves, in- 
dicating this is a starburst galaxy. CDFS-6482 has a rel- 
atively high [OIII] / H/3 value (lower limit), but [Nil] /Ha 
is not so high. This high [OIII] / H/3 ratio may not be 
due to the presence of AGN but due to a higher gas 
temperature of the HII regions, as has been p ointed out 
for som e high-redshift star-forming galaxies (|Erb et al.l 
l2006a| ). Indeed, there are few local AGNs distributed 
in the range of log([NI I]/Ha) < —0.5 (e.g., Figure 3 in 
IFinkelstein et al.ll201lD. 

IFinkelstein et al I (120091) have also classified a similar 
object to ours as a star-forming galaxy (blue triangle in 
Figure [3]). A similar argument is made for COSMOS- 
13636 and COSMOS-30679 whose [Nil] /Ha ratios are 
also modest. In contrast, COSMOS-43982 has a very 
high [Nil] /Ha ratio. Because we infer that this object is 
more likely to be an AGN than a star-forming galaxy, we 
do not use this object in the following discussion. 



1.5 




-1.5 -1 -0.5 
Log([NII]/Ha) 



Fig. 3. — BPT diagram. The solid curve represents the 
boundary between star-forming galaxies and AGNs proposed by 
IKewlev et "all if^OOlD . while the dotted curve is the boundar y em- 
pirically defined using SDSS objects JRauffmann et al. 2003]). The 
red circles denote our LAEs. For the purpose of display, COSMOS- 
13636 and COSMOS-43982, whose [OIII] / H/3 is not constrained, 
are placed at log([OlII] / H/3) = 0.0. The red triangles, blue cir- 
cles, and blue tr i angle show, respectively, two z ~ 2. 3 LAEs in 
IFinkelstein et al.1 <201lD . four bins of z ~ 2.2 LBGs inlErb et all 
l|2006al ). and a lensed LBG at 2 = 2.73 IjFinkelstein et al.N2009T l 



3. SED FITTING 

We perform SED fitting to our objects to derive dust 
extinction values, SFRs, and stellar masses. The proce- 
dure of the SED fitting is the same as that of lOno et al.l 
(|2010al) . For the CDFS objects, we use 12 bandpasses: 
B, V, R, I, z, J, H, K data from the MUSYC public data 



release E3 (|Cardamone et al.|[2010t) . and Spitzer/IRAC 
3.6, 4.5, 5.8, and 8.0 /im photometry from the Spitzer 
legacy survey of the UDS field. For the COSMOS ob- 
jects, we use 11 bandpasses: B,V,r',i', and z' data 
taken with Subaru/Suprime-Cam, J data taken with 
UKIRT / WFC AM K s d ata taken with CFHT / WIRC AM 
(jMccracken et al J 12010) . and Spitzer/IRAC 3.6, 4.5, 5.8, 
and 8.0 /jm photometry from the Spitzer legacy survey 
of the UDS field. We use neither u*/U nor NB387-band 
data, since the photometry of these bands is contami- 
nated by the IGM absorption and/or Lya emission. Ta- 
bles [3] and 0] summarize the broadband photometry of our 
objects. Basically, the uncertainties in the optical pho- 
tometry include both photometric errors and systematic 
errors associated with aperture correction and photomet- 
ric calibration. 

We use t he stellar population sy nthesis model of 
GALAXEV (IBruzual fc Chariot] 12003) in cluding nebu- 
lar emission (jSchaerer fc de Barrosl 20091). an d adopt a 
Salpeter initial mass function (lSafpeterlll955l) . Because 
LAEs are metal poor star-forming galaxies, we choose 
constant star formation models w ith a metallicity o f 
Z/Z Q = 0.2. We use Calzetti's law (jCalzetti et alJl200O) 
for the stellar con tinuum extinction, E(B — V)„. We ap- 
ply iMadaul (|1995| ) 's prescription to correct for the IGM 
attenuation; at z ~ 2.2, continuum photons shortward of 
Lya are absorbed by 18 %. Figure [4] shows the best-fit 
model spectra with the observed flux densities for in- 
dividual objects. The physical quantities derived from 
SED fitting arc summarized in Table [5] Because the 
continuum emission of COSMOS-30679 is blended with 
a foreground obje ct, we deblend the source wi th GALFIT 
(|Peng et al .1120021: see lNakaiima et al.|[2012bl for details). 
While photometry data both before and after deblending 
are given in Table |4j we use the latter for SED fitting. 

4. PHYSICAL QUANTITIES OF OUR LAES 

4.1. Velocity Offset Between Lya and Nebular Lines 

We calculate a quantity called the velocity offset of the 
Lya line: 



•^Lya 



(3) 



where ZLya is the redshift of Lya emission and z sys is the 
systemic redshift given in Table [2j 

We show the Lya and Ha profiles in Figure [5] Most 
objects have a Lya profile which is asym metric and/or 
double peaked. Similar to previous studies (jSteidel et al.1 
[2^lltlYaii7et^[201ll) , we define -^Lya 

to be the redshift 

corresponding to the wavelength at the highest peak. 
The peak wavelength is determined as follows. First, we 
roughly constrain the peak position by the eye. For the 
CDFS objects, we cannot identify where is the highest 
peak, since the second peak is within the la error in the 
height of the first peak. In these cases, we regard both 
as peak candidates and include their wavelength differ- 
ence as the error in z-^ ya (as listed in Table [T]). Then, 
we fit a Gaussian to the profile only around the peak to 
derive the peak wavelength, in order to avoid systematic 
effects due to the asymmetric profile. The Lya line of 
COSMOS-30679 is partly contaminated by a cosmic ray. 
It is, however, unlikely that there is a true peak under 

15 http://www.astro.yale.edu/MUSYC/ 
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TABLE 3 

Broadband Photometry of our Sample (CDFS) 



Object 


B 


V 


R 


I 


z 


J 


H 


K 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


CDFS-3865 


23.01 
(0.01) 


22.94 
(0.01) 


22.92 
(0.01) 


23.14 
(0.07) 


22.93 
(0.09) 


22.73 
(0.18) 


22.27 
(0.12) 


22.38 
(0.23) 


22.82 
(0.05) 


22.82 
(0.08) 


22.51 
(0.32) 


23.00 
(0.56) 


CDFS-6482 


23.93 
(0.02) 


23.87 
(0.03) 


23.78 
(0.03) 


23.95 
(0.14) 


23.67 
(0.16) 


23.50 
(0.34) 


23.36 
(0.31) 


23.07 
(0.39) 


22.88 
(0.05) 


22.83 
(0.08) 


23.34 
(0.60) 


99.99 
(-) 



Note. — All magnitudes are total magnitudes. 99.99 mag indicates a negative flux density. Magnitudes in 
parentheses are lcr uncertainties. 



TABLE 4 

Broadband Photometry of our Sample (COSMOS) 



Object 


B 


V 


r> 


%> 


z' 


J 


K a 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


COSMOS-13636 


24.43 


24.21 


24.35 


24.19 


24.24 


23.10 


23.43 


24.10 


23.75 


99.99 


99.99 




(0.01) 


(0.03) 


(0.03) 


(0.04) 


(0.09) 


(0.34) 


(0.26) 


(0.38) 


(0.53) 


(-) 


(-) 


COSMOS-30679" 


24.05 


23.12 


22.91 


22.46 


22.33 


21.15 


21.83 


22.12 


22.57 


99.99 


23.06 




(0.01) 


(0.01) 


(0.01) 


(0.01) 


(0.02) 


(0.07) 


(0.07) 


(0.07) 


(0.21) 


(-) 


(2.54) 


COSMOS-30679 6 


24.76 


23.82 


24.44 


24.09 


23.49 


22.31 


23.29 












(0.03) 


(0.11) 


(0.28) 


(0.30) 


(0.27) 


(0.27) 


(0.28) 


(-) 


(-) 


(-) 


(-) 



Note. — All magnitudes are total magnitudes. 99.99 
parentheses are 1<t uncertainties. 
a Before deblcnding of a foreground source. 
b After deblcnding of a foreground source. 

TABLE 5 
Results of SED fitting 



mag indicates a negative flux density. Magnitudes in 



Object 
(1) 



E(B-V)* 
(2) 



SFRsed 

(Mq yr" 1 ) 
(3) 



(10 9 Mq) 

(4) 



(5) 



CDFS-3865 
CDFS-6482 
COSMOS-13636 
COSMOS-30679 



1-0.009 
-0.009 
H0.026 
-0.018 
,+0.018 
-0.079 
u.ozo_ 026 



0.185] 
0.185j 
0.273j 



312 
83 
131F 
7510 4 



+ 41 
-43 
H34 
-18 
-17574 
-1244 
201248 
3097 



3.18 
5.30 
1.99 
19.75 



+0.21 
0.13 
+ 1.18 
0.80 
+ 0.39 
1.06 
-6.53 
-5.80 



17 

5 

21 
20 



Notes. Stellar metallicity is fixed to 0.2 Zq. 

(1) Object ID; (2) Dust extinction; (3) Star formation rate; (4) 
Stellar mass; (5) x 2 °f t nc fitting. 

the cosmic ray. Even if there is a flux peak just under 
the cosmic ray, our discussion remains unchanged. 

We obtain Av Lya = 190t?| km s" 1 (CDFS-3865), 
75±H km s- 1 (CDFS-6482), 18 ± 16 km s" 1 (COSMOS- 
13636), and 170 ± 16 km s" 1 (COSMOS-30679). Thus, 
three out of the four have a positive Avi,y a beyond the 
2a uncertainty. 

There are four LAEs in the literature which have both 
a Lya spectrum and a systemi c-velocity measu r ement 
from a nebula r line: two from iMcLinden et al.1 (|2011l ) 
and two from iFink clstein e t al.l (|2011D . The AuL yQ of 
these four LAEs ranges from ~ 100 to ~ 300 km s _1 , 
which is similar to our Av-^ ya measurements. We com- 
bine these four LAEs with our four to construct a sam- 
ple of eight LAEs at z ~ 2 — 3, and investigate av- 
erage physical properties of LAEs using them. Fig- 
ure [5] presents the histogram of AwL yQ for the eight 
LAEs. All eight have a positive but relatively small 
A^Lya of up to ~ 300 km s _1 . The average of the eight 
is A^Lya = 175 ± 35 km s _1 , which is systematically 
smaller, by ~ 200 - 300 km s -1 , than that of LBGs, 



Av h y a ~ 400 km s" 1 (IPettini et all 12002 IShaplev et al.1 
[20031 ISteidel et al] 120101: IRakic et all 120111 ). at similar 
redshifts of z ~ 2-3. There ar e also two Lya blobs 
(LABs; e.g.. ISteidel et all [20001: iMatsuda et all 120041: 
lYang et a l. 2009) at z ~ 2.3 whose Av LyQ is measured 
(jYang et al.l 120 111) . Both have Au LyQ = - 200 km s" 1 
which is comparable with our measurements for LAEs 
but smaller than those for LBGs. 

Figure [7] shows the rest-frame Ly a EW as a func- 
tion of AuLva, wh ere data for LBGs dReddv et al]|2008t 
ISteidel et all2010l) and LABs (|Yang et al.ll201lD are also 
included in order to cover a wide baseline of EW, — 200 
A. We fit a linear function to the data points (dotted line 
in Figure [7]), and find an anti-correlation between Avi^ ya 
and Lya EW that EW(Lya) decreases with increasing 
Ai>L y a- It is known that there is an anti-correlation be- 
tween AfLyq - Af a b s and EW(Lya) in LBGs at z ~ 3 
(jShaplev et al.|[2003l) . However, it should be noted that 
our finding is different from that, because we measure 
the velocity offset of Lya directly from the systemic ve- 
locity, while At>L yQ ; - Aw a bs depends on the kinematics of 
gas outflows and the Lya emission mechanism. 

4.2. Velocity Offset Between LIS Absorption Lines And 
Nebular Lines 

We examine LIS absorption lines using Magcl- 
lan/MagE echcllette spectrograph data with a medium- 
high spectral resolution of R ~ 4100 (M. Rauch et al., in 
preparation) . We inspect one order of the echellc spectra 
ranging from 3788 A to 4419 A in the observed frame. 
This wavelength range covers Si II A1260, O I A1302, Si 
II A1304, and C II A1335 lines. These LIS absorption 
lines are associated with the neutral interstell ar medium 
(ISM: IPettini et al.ll200l IShaplev et an[2003l) . 

Since none of these lines is identified in the individ- 
ual spectra because of the too faint continua, wc make 
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Fig. 4. — Upper (lower) panel shows the SEDs of CDFS (COSMOS) objects. The filled squares denote the photometry points used for 
SED fitting, while the open squares are those omitted in SED fitting due to the contamination of Lya emission and IGM absorption. The 
red lines present the best-fit model spectra, and the red crosses correspond to the flux densities at individual passbands expected from the 
best-fit models. 



a composite spectrum of the four objects as follows. For 
each object, we shift individual spectral data from the 
observed to the rest frame using z sys given in Table [2j 
Then, we stack the spectra of the four objects with sta- 
tistical weights based on the S/N ratios at 1250-1340 
A. 

Figure |S] presents the composite spectrum. An absorp- 
tion feature is seen near each of the four lines' wave- 
lengths. For each feature, we fit a Gaussian profile to 
obtain the rest-frame equivalent width, EW(LIS), and 
the velocity offset of the line center with respect to the 
systemic velocity, Au a b s , as summarized in Tabled We 
securely detect the Si II A 1260 line and the blended O I 
A1302 + Si II A1304 lines at the 5.2 a and 3.7 a levels, 
respectively, and marginally detect the C II A1335 line at 
the 2. Oct level. Note that all lines are blue-shifted with 
respect to the systemic velocity. 

The weighted mean offset velocity of these absorp- 
tion lines is Av a b s = —179 ± 73 km s _1 . This value 
is comparable with th ose of LBGs, which are typically 

Av ah s 150 km s" 1 (jShaplev et al.H2003t ISteidel et all 

120101) . Thus, LAEs and LBGs have similar Av a b s values, 
in contrast to the significant difference in AuL yQ . 

4.3. Dust Extinction 



TABLE 6 
Interstellar Absorption Features 



Ion A ros t / EW(LIS) cr Au abs 

(A) (A) (A) (kms" 1 ) 
(1) (2) (3) (4) (5) (6) 

Si II 1260.4221 1.007 -1.21 0.23 -162 ±95 

Oi 1302.1685 0.04887 -1.02< 7 ) 0.28( 7 > -174±180( 7 > 

Si II 1304.3702 0.094 -1.02< 7 ) 0.28( 7 > -174±180( 7 > 

C II 1334.5323 0.1278 -0.52 0.25 -209 ±110 

(1) Absorption line; (2) Rest-fra me vacuum wavelength; (3) Tran- 
sition oscillator strength (sec e.g.. Pettini et al. 2002; Shaplcv et al. 
2003); (4)-(5) Rest-frame EW and its 1 cr; (6) Velocity offset. (7) 
Values for the blended O I A1302 and Si II A1304 lines assuming 
that the central rest-wavelength wavelength is A = 1303.2694 A. 



The color excess of the stellar continuum obtained 
from the SED fitting is E(B - V)* = 0.185±g;gg| 
(CDFS-3865), 0.185±g;g?| (CDFS-6482), 0.273±g;{£§ 
(C OSMOS-13636) and 0.528±g;{gg (COSMOS-30679) 
for iCalzetti et all (|2000D 's extinction law. These val- 
ues are comparable to those of relatively faint (stacked) 
LAEs at z ~ 2, E(B - V)* = 0.27^;^ (|Nakaiima et all 

l2012ah and E(B - V)* = 0.22t°;™ (|Guaita et alll20Tl . 
but are larger than those of the two bright LAEs observed 
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Fig. 5. — Comparison between Lya (blue line) and Ha (red line) emission profiles. The solid and dotted lines indicate the systemic 
and Lya redshifts, respectively. Top left: CDFS-3865 with Ai> Lya = 190lfg km s _1 . Top right: CDFS-6482 with Av Lya = 75±H km s _1 . 
Bottom left: COSMOS-13636 with Av Lya = 18 ± 16 km s _1 . Bottom right: COSMOS-30679 with Av hya = 170 ± 16 km s" 1 . 

Since the spectra of CDFS-3865 and CDFS-6482 cover 
H/3 as well as Ha, EE we measure the Balmer decre- 
ment (Ha/H/3) to constrain the color excess of nebu- 
lar emission, E(B — V) &as . The decrement values ob- 
tained are 2.96 ± 0.87 for CDFS-3865 and > 1.74 (using 
H/3's 2a limit) for CDFS-6482, which are converted into 
E(B - y) gas = 0.03tg;|g and E(B - T/) gas > 0, respec- 
tively , ad opting the intr insic ratio of 2.86 (|Osterbrockl 
Il989h and lCalzetti et all (|2000l )'s extinction law. 

Figure M plots E(B - V)* vs. E(B - V) gas for these 
objects. CDFS-3865 has E(B - V) gas ~ E(B - V),. 
Similarly, CDFS-6482's lower limit of E(B - V) gas is con- 
sistent with its E(B — V)*. Thus, hereafter, we assume 
E(B - V) = E(B - V0 gas ^ E(B - V)* as proposed by 
lErb et all (|2006b| ) for starburst galaxies. 

4.4. Star Formation Rate 

The Ha luminosity is thought to be a reliable indicator 
of SFR. We calculate the SFRs of o ur LAEs from the ir 
dust-corrected Ha luminosities using iKcnnicuttl (|1998l )'s 
formula: 

SFR(M Q yr- 1 ) = 7.9 x 10" 42 L(Ha) (erg s" 1 ), (4) 

where a Salpeter IMF is assumed. We obtain SFR = 
19O±13Af yr" 1 (CDFS-3865), 48t^°M Q yr" 1 (CDFS- 
6482), 17±^M yr" 1 (COSMOS-13636), and 45 ± 
5M Q yr" 1 (COS MOS-30679). The s e valu es are compa- 
rable to those of iFinkelstein et aTl (]201lD 's two LAEs, 

16 COSMOS-30679 also has H/3 data, but the line is heavily 
contaminated by a sky emission line. 
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Fig. 6. — Histograms of AiiLya for the eight LAEs studied in 
this paper (red) and 41 LBGs given by Stcidel et al. (2010) (blue). 



bv IFinkelstein et alJlpOll . E(B-V)* = 0.09±0.05and 
E (B — V)* = 0.10 ±0.10. Since our objects are as bright 
as IFinkelstein et all (|2011l )'s. this difference may indicate 
that LAEs at z ~ 2 have a wide range of dust extinction 
even at the same Lya luminosity range. 
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Fig. 9.— E(B - V) g & B plotted against E(B - V)* for oui two 
LAEs with a Balmer decrement measurement. The d otted and 
dashed lines correspond to E(B — V) t = E( B — V)„ ag l|Erb et alj 
|2006rj ) and E(B - V), = 0AAE(B - V) gas l lCalzetti et aI.II2000n , 
respectively. 



Fig. 7. — Rest-frame EW(Lyo) plotted against Av^ ya ,. The 
red circles are our LAEs. The red squar e and the red triangle 
show the LAEs by McLinden et al. (2011) and Finkclstcin et al. 
I|2011l ). respectively. The blue symbol indicates the average of 41 
LBGs, with the error b ars corresponding t o the 68 percentiles of 
the Ai)L yc « distribution (Stcidel et al. 2010) and the EW distribu- 
tion l|R cddv ct al. 2008]). The magenta circles denote the LABs by 
lYang et al.l H2011I ). The dotted line is the best-fit linear function 
to all the data points. 
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Fig. 8. — Composite FUV spectrum of the four LAEs. The 
spectrum has been normalized to unity in the continuum. The 
dotted lines indicate the rest-frame vacuum wavelengths of four 
LIS absorption lines. The spectrum is plotted with 0.5 A pix~ 1 
sampling so that the profile of the blended O I and Si II lines be 
clearly shown. 

36.8 ± 5.8 M Q yr- 1 and 23.6 ± 6.5 M G yr' 1 , derived 
from the Ha luminosity, but are larger than the average 
values of z ~ 2.2 L AEs obtained from stacking analysis, 
10 - 2 M Q yr" 1 (jNilsson et al.l 120111: iNakaiima et al.l 
l2012al) . This difference is reasonable, because LAEs with 
Ha detection are generally brighter than average LAEs. 
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Fig. 10. — SFRs from SED fitting and from the Ha luminosity. 
From left to right, plotted are COSMOS-13636, COSMOS-30679, 
CDFS-6482, and CDFS-3865. 

Figure \W\ plots the SFR from SED fitting against 
that from the Ha luminosity for our objects. Although 
the two SFR values are comparable for CDFS-3865 and 
CDFS-6482, they are largely different for COSMOS- 
13636 and COSMOS-30679, indicating that models with 
very different star formation histories can fit the observed 
SED almost equally well. This figure thus demonstrates 
that it is important to measure the Ha luminosity to de- 
rive a reliable SFR. In the rest of this paper, we quote 
the SFRs obtained from the Ha luminosity. 

4.5. Lya Escape Fraction 

The Lya escape fraction of a galaxy, f^ a , is defined 
as the ratio of the observed Lya flux to the intrinsic Lya 
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flux produced in the galaxy. This quantity can be a probe 
of the distribution and kinematics of the ISM. For exam- 
ple, the outflowing ISM would make larger since 
the number of resonant scattering and thu s the chance 
of absorption by dust are reduced (e.g., iKunth et al.1 
119981 lAtek et all l2()()8li. A clumpy distr ibution of the 
ISM would also make fiZ a larger (e.g., lHansen STOhl 
120061 iFinkelstein etaD f2008h . In the clumpy ISM, HII 
regions where both Lya and nebular lines originate are 
surrounded by clumpy gas cl ouds with well- mixed HI gas 
and dust (see, e.g., Fig. 1 of lNeufeldlH990l ). We discuss 
the extent of dumpiness for our objects in M5.3.21 

Assuming the Case B recombination whe re the intrin- 
sic Lya/Ha ratio is 8.7 ( Brocklehurst1 ll97lh . we estimate 
the Lya escape fraction as: 



Lya _ 
esc — 



Lobs (Lya) _ L obs (Lya) 



L;„t (Lya) 8.7L int (Ha) 



(5) 



where subscripts 'int' and 'obs' refer to intrinsic and ob- 
served quantities, respectively, and Li nt (Ha) is obtained 
by correcting the observed Ha luminosity for dust ex- 
tinction. We derive the L ODS (Lya) of each object from 
its narrow and broad-band photometry on the assump- 
tion of a flat continuum, taking accout of the exact posi- 
tion of the Lya emission in the NB387 response function 
and applying Madau (1995) 's prescription to correct for 
the IGM attenuation. Note that these f^. a values are 
not measurements but estimates, because we assume the 
Case B recombination. 

We find f^/ c a = 0.14 ± 0.03 (CDFS-3865), 0.29 ± 0.20 
(CDFS-6482), 0.57±0.13 (COSMOS-13636), and 0.17 ± 
0.03 (COSMOS-30679). These values are much higher 
than t he average value of z ~ 2 star-forming galaxies, 
~ 5% (jHaves et al.ll2010D . but similar to the median value 
of 89 L AEs at 2 < z < 4, ~ 29%, obtained bv lBlanc etlll 
()2011|) . It is interesting because our LAEs have relatively 
large E(B-V) values (Tabled . This would suggest that 
some mechanisms allow Lya photons to escape from the 
moderately dusty ISM (jAtek et alJl2008T ). 

4.6. Size 

Assuming that stars dominate the total mass of the 
luminous (i.e., Ha emitting) part of galaxies, we infer 
the size of our objects from the virial theorem as: 



G 



AL 



M*/10 10 M Q 



CTv(Ha) 2 



(cr v (Ha)/100 km s" 1 ) 2 



(6) 

where M* is the stellar mass and <7 v (Ha) is the velocity 
dispersion measured from the Ha line, corrected for the 
instrumental velocity dispersion (<7i ns t = 86 [91] km s _1 
for NIRSPEC [MMIRS]). It is known that for local galax- 
ies stellar masses dominat e total masses within effe ctive 
radii (see e.g., Fig 2. a in Ivan Dokkum et al.l I2009f ) and 
that for high-z galaxies total masses c an be larger than 
stellar masses up to a fac tor of several (|Erb et alj|2006bt 
van Dokkum et ail 120091 ) . Thus, our size estimates ap- 
pear to be reasonable approximations of the true values 
within a factor of two or so. 

We obtain r = 1.1±H kpc (CDFS-3865), 5.0tJ;£ 
kpc (CDFS-6482), 9.0±^| kpc (COSMOS-13636), and 
107.5±|f:| kpc (COSMOS-30679). Similarly, we derive 
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Fig. 11. — Size plotted against stellar mass. The red circles de- 
note our LAEs, while the blue circles represent 20 LBGs which have 
reliable measurements of stellar mas s and He? velocity dispersion 
HErb et al.M2006b1 ; ISteidel et aHl20T0D . 

sizes for 20 LBGs of lErb et al.l (|2006bf ) which have mea- 
surements of stellar mass and Ha velocity dispersion. 
Here we multiply the stellar masses of Erb et al.'s ob- 
jects by 1.8, since they assume a Chabrier IMF (jChabrierl 
120031) . As found from Figure [TT] the LAEs and LBGs 
roughly lie on a single size-stellar mass relation, sug- 
gesting that typical LAEs are smaller than LBGs just 
because they are less massive. 

The physical quantities presented in this section are 
summarized in Table [7] 

5. DISCUSSION 
5.1. A VLya and the Gas Motions of LAEs 

Due to the resonant nature of Lya, the observed Lya 
line of a galaxy has a complicated profile depending on 
the kinematics and geometry of the ISM. A Lya source in 
a simple static gas cloud produces a symmetric double- 
peaked profile centered at 1216 A due to significant res- 
onan t scattering at 1216 A (|Harringtonl 119731 : iNeufeldl 
119901: IDiikstra et ail 120061) . If the bluer peak is heavily 
absorbed by the intervening IGM along the line of sight, 
only the redder peak will be observed. In the case of 
outflow, the Lya emission line should s how an asymmet- 
ric pro file similar to a P Cygni profile. iVerhamme et all 
(|2006| ) have explained some of observed Lya profiles with 
a strong red peak and a weak blue peak with mod- 
els of an expanding shell that absorbs Lya photons at 
around 1216 A, although the surface brightness distribu- 
tion may not be explaine d by such wind shells (see, e.g., 
IBarnes fc Haehneltll2009D . In short, when an outflow ex- 
ists, the observed Lya line should show an asymmetric 
profile with a strong peak rcdshiftcd with respect to the 
systemic velocity. 

We measure the offset velocity of the Lya peak, 
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A^Lya, in N4.ll To quantify the asymmetry of the 
Lya line of our objects, w e use the weighted skewness, 
S w ([Kashikaw a et all 12006ft . Skewness is the third mo- 
ment of the flux distribution, and the weighted skew- 
ness is defined as the product of the skewness and the 
line width. A positive S w means that the Lya pro- 
file has a red tail as in the case of outflow. We obtain 
S w = 6.03±0.51 (CDFS-3865), 5.45±1.59 (CDFS-6482), 
1.01±0.76 (COSMOS-13636), and 5.31±4.06 (COSMOS- 
30679) EH 

As found in ij4.1[ three out of the four objects have a 
positive A^Lya beyond the 2a uncertainty. This result, 
combined with the positive measurements of S w , leads us 
to a conclusion that the Lya emission of our objects is 
mostly originated from outflowing gas. This conclusion 
is supported by the fact obtained in M4.2l that the LIS ab- 
sorption lines of the composite spectrum are blue-shifted 
with respect to the systemic velocity. 

Using sim ulations of Lya radiative trans fer in outflow- 
ing galaxies, IVerhamme et a l. (2006, 2008]) suggest that 
not only LIS absorption lines but also the Lya line can 
be used to estimate the outflow velocity. Their simula- 
tions assume that a galaxy is surrounded by a spheri- 
cally symmetric shell-like outflowing gas cloud in which 
Hi gas and dust are well mixed. They find that for rel- 
atively low Hi column densities of Ah < 10 20 cm -2 , the 
peak of the Lya profile emerges near the outflow velocity, 
giving AvLya ~ v out , while for high column densities of 
Ah > 10 20 cm~ 2 , the peak is offset twice the outflow ve- 
locity, giving AwLya ~ 2u out . In any case, it is suggested 
that A«Lya positively correlates with v on t- 

On a reasonable assumption that A«Lya positively cor- 
relates with the speed of an outflow within a factor of 2, 
one can test the hypothesis that Lya photons can escape 
more easily for a larger outflow velocity, d ue to the re- 
duced number of resonant scattering (e.g., iKunth et al.l 
[1991) . If this is the case, we should find a positive cor- 
relation between AwL yQ and EW(Lya). However, Figure 
[7] indicates an opposite tendency. Furthermore, we find 
]that LAEs and LBGs have similar average Au a b s , 



175 km s 1 , these studie s may be overestimatin g a-jji- 



i.e., similar average outflow velocities. These two findings 
would suggest that outflows are not the physical origin 
of large-EW(Lya) objects. 

In ij4.ll we find that the average AuL yQ of LAEs is 
175 ± 35 km s _1 , which is significantly smaller than that 
of LBGs (~ 400 km s _1 ) at similar redshifts. This re- 
sult is important not only for understanding the physical 
origin of Lya emission in galaxies, but also for prob- 
ing cosmic reionization with LAEs. If LAEs at z > 6 
have similarly small AuL yct values, the amount of Lya 
photons scattered by the IGM, as used to constrain the 
epoch of reionization, may be in need of revision. For 
example, ISantosI (j2004j ) has examined the transmission 
through the IGM of Lya photons emitted from a galaxy 
for two cases, AvLya = and 360 km s _1 , the latter of 
which is comparable to the average AwLya of z ~ 2 - 
3 LBGs. Some recent reionization studies using LAEs 



fraction of the IGM, imi at z > 6 (e.g., 


Kashikawa et al. 


2006 


:lOta et al.ll2008HOuchi et aU201d 


Kashikawa et al. 


2011 


). However, if z > 6 LAEs have AwL yQ as small as 



17 The S w of COSMOS-30679 is calculated after masking the 
wavelength range affected by the cosmic ray. 



Similarly. lOno et al.l ([20121) . iSchenker et all poll , and 
iPentericci et al.l ( 2011 ) have derived xm as large as ~ 40 
- 60% from a significant drop in the fraction of large- 
EW(Lya) galaxies from z ~ 6 to 7. If such a high value 
were correct, reionization would take place very late, 
which cannot easily be reconciled with constraints fr om 
the Lyman alpha forest opacity (jBecker et al.l 120071) or 
the large value of the Thomson op tical depth, t = 0.09 , 
obtained by WMAP observations (jDunklev et al.l 120091: 
iKomatsu et al.ll201lD . Future Lya emission models for 
reionization studies would need to use our result of a 
small average AfLy Q value if A«Lya does not significantly 
evolve over cosmic time. 

5.2. Correlations Between Outflow Velocity and 
Physical Properties 

5.2.1. Outflow Velocity and SFR 

We first examine the correlation between A«Lya, which 
represents the outflow velocity within a factor of two as 
discussed in «j5. It and the SFR for LAEs and LBGs. For 
the present-day universe, iMartinl ([2005t ) has found m a 
sample of ULIRGs, LIRGs, and starburst dwarfs that 
those with a larger outflow velocity tend to have a higher 
SFR, roughly following a power- law of u ut oc SFR 35 . 
Note that this sample spans four orders of magnitude 
in the SFR. This correlation implies that galaxies with a 
high SFR tend to have more massive stars and SNe which 
drive outflows. As pointed out by some authors (e.g., 
lMartin|[2005l ; iSteidel et al.ll2010D . the correlation flattens 
and is difficult to see for objects with intermediate SFRs 
of 10-100 M Q yr^ 1 . 

From the left panel of Figure Q21 there may exist a posi- 
tive correlation between AwLya and SFR for LAEs, which 
is roughly consistent with A^Lya oc SFR, although the 
statistics is not ver y good. For LBGs, which are take n 
from the literature ((Erb et al.ll2006bHSteidel et al.ll2010( ). 
a much flatter correlation is seen. Thus, it is found 
that LAEs and LBGs do not follow the same power 
law, and that LAEs have comparable to or systemat- 
ically smaller outflow velocities at a given SFR, even 
when we take account of the possibility that A^Lya may 
overestimate w ou t up to a factor two. Theoretically, the 
slope of the correlati on depends on the dominant phys- 
ical outflow process (jKornei et alJ[2012T ). A linear cor- 
relation, v out oc SF R, is expected for r adiation pressure 
dominant outflows (Sh arma et al.ll201l|) . while shallower 
slopes are indicative o f ram pressure dominant outflows 
(jHeckman et al.ll2000r ). It appears that outflows in LAEs 
are at least consistent with the radiation pressure case. 

5.2.2. Outflow Velocity and a v (Ha) 

We then examine how AwL ya correlates with er v (Ha), 
whi ch is rel a ted to the gravitational potential of galax- 
ies. IMartinl (|2005l ) has found that the terminal outflow 
velocity always approaches the galactic escape velocity, 
which is consistent with a theoretical prediction (e.g., 
iMurrav et all I2005D . The middle panel of Figure [12] 
presents AvLya as a function of ay (H a) for our LAEs , 
together with tho se of ~ 20 LBGs ([Erb et al.l l2006bt 
ISteidel et al.l 120 lOh . It is not clear whether there exists 
a correlation with the given poor statistics and large er- 
rors. However, LAEs are found to have comparable to or 
smaller outflow velocities than LBGs at a given er v (Ha). 
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5.2.3. Outflow Velocity and SFR Surface Density 

Finally, we examine the correlation with the SFR sur- 
face density, Esfr, which represents the intensity of star 
forming activity. The right panel of Figure [12] plots 
A^Lya against Esfr for our LAEs and LBGs taken 
from the literature. Here we estimate Esfr by divid- 
ing the SFR by irr 2 , where r is the size derived in fl4.6l 
The Esfr values of our objects are 157t|i yr -1 
kpc- 2 (CDFS-3865), 1.9±gf M Q yr" 1 kpc" 2 (CDFS- 
6482), 0.22lg;S 8 Me yr -1 kpc" 2 (COSMOS-13636), and 
3.89± 2 j 2 x 10" 3 M yr" 1 kpc" 2 (COSMOS-30679). 
Again, LAEs are found to have comparable to or smaller 
outflow velocities than LBGs at a given Esfr. 

To summarize, we find that LAEs tend to have compa- 
rable to or smaller outflow velocities than LBGs and that 
LAEs and LBGs have different slopes of the At>Ly Q -SFR 
relation. These findings imply that the physical origin 
of LAEs' outflows appears to be different from that of 
LBGs'. 

5.3. Why LAEs Have Strong Lya Emission? 

In this subsection, we discuss the physical origin of 
strong Lya emission from LAEs. There are various mech- 
anisms capable of enhancing Lya emission: shock heat- 
ing by an inflow or outflow, very weak dust extinction, a 
peculiar geometry of dust /gas clouds, a low neutral hy- 
drogen column density, and a high outflow velocity. The 
high outflow velocity hypothesis is found to be unlikely 
in ^5.11 (see also the left panels of Figure [T5|) . Here we 
examine the weak-dust extinction hypothesis, clumpy- 
cloud hypothesis, and low- Ah hypothesis. 

5.3.1. Weak Dust Extinction 

It is possible that Lya photons survive through reso- 
nant scattering in gas clouds with little dust. In this case, 
one should find significantly small E(B — V) in LAEs and 
a negative correlation between E(B — V) and EW(Lya). 
The middle panels of Figure [TJ plot EW(Lya) and f^g a 
as functions of E(B — V). There is no s igni ficant correla- 
tion in either panel. lHaves et al.l (I2010D and lKornei et al.l 
(|2010t ) have reported an anti-correlation between f^/ c a 
and E(B — V) in Ha emitters at z ~ 2.2 and LBGs at 
z ~ 3, respectively. However, in our f^ a vs E{B — V) 
plot, there are objects with relatively high f^/ c a of 0.3-0.6 
even in a moderately high extinction of E{B — V)= 0.2- 
0.3. Therefore, it is not clear if weak dust extinction 
solely can explain the strong Lya emission. 

5.3.2. Clumpy Clouds 

The gas distribution in LAEs may not be smooth and 
spherically symmetric, but clumpy. In a clumpy geome- 
try, dust grains are shielded by HI gas, and Lya photons 
are resonantly scattered on the surfaces of clouds with- 
out b eing absorbed by dust (|Neufeldlll99ll lH anscn &~Ohl 
120061 ). Because continuum photons are absorbed through 
dusty gas clouds, the ratio of Lya to UV continuum 

18 Although many studies have adopted Petrosian radius for the 
size of a galaxy, it is known that Petrosian radius is likely to over- 
estimate the area of the star forming region of galax ies especially 
for hi gh- z objects with a clumpy star forming region l IKornei et al.l 



fluxes, or EW(Lya), is enhanced. For a further discus- 
sion of the clumpy-cloud hypothesis, we calculate the 
dump iness parameter, q, introduced bvlFinkclstci net al.l 
((200l : 

Q = 7-(Lya)/Ti2i 6 , (7) 

where r(Lya) and 71216 arc defined as e~ T ( Lya ^ = 
L obs (Lya)/^ (Lya) and e" 7 " 1216 = 10"° -^i 216 e{b-v) 
with t he extinction coeffic ient at A = 1216 A, fci2i6 = 
11.98 IjCalzetti et al.l 12000? ). This par ameter is used to 
diagnose the geometry of the ISM (e.g..lFinkelstein et al.l 
I2008t iKornei et al.ll2010l : iNakaiima et al.ll2012al) . If the 
geometry of the ISM is clumpy, q is smaller than unity. 
In the case of q > 1, Lya photons are more preferen- 
tially absorbed by dust through the relatively homoge- 
neous ISM. 

The q values of our objects are calculated to be 
q = 0.96 ± 0.11 (CDFS-3865), 0.61 ± 0.34 (CDFS-6482), 

0. 19±g;g£ (COSMSO-13636), and 0.30±°;° 3 (COSMOS- 
30679); three out of the four clearly have q < 1 beyond 
th e uncertainties. These r esults are consistent with that 
of INakaiima et all (|2012aD who obtained q = 0.7±°[\ for 
z ~ 2.2 stacked LAEs. On the other hand, some stud- 
ies have found more LAEs with q > 1 at z ~ 2 (e.g., 
iNilsson et all [20091 : [Haves et al.|[2010t) . 

The right panels of Figure [TBI show and EW(Lya) 
as functions of q for the LAE sample. There may exist 
a negative correlation between q and f^. a , although not 
very clear due to the poor statistics. 

5.3.3. Low N H 

Finally, we examine the possibility that the Ah of 
LAEs is relatively low. To constrain Ah for LAEs and 
L BGs, we use the dependence of Ai>L yQ on Ah found 
bv lVerhamme et al.l (|2006l . I2008f ). As mentioned in fUj 
they find A« LyQ ~ 2w out for A H > 10 20 cm -2 while 
A^Lya ~ v out for Ah < 10 20 cm -2 ; Av^ ya is smaller for 
lower Ah due to a smaller amount of wavelength shift by 
resonant scattering. Since |Aw a b s | is equal to v ou t in an 
expanding shell, one can thus use AwL ya and Ai> a b s to 
distinguish between high and low Hi column densities. 

The LAEs we study have average velocities of Av^ya = 
175 ± 35 km s" 1 and Aw abs = -179 ± 73 km s" 1 , 

1. e., A^Lya ~ |Ai? a b s |, while LBGs have At>L V q ~ 2- 
3 x |Aw abs | (jSteidel et al.ll20Tol : iRakic et al.|[201l[) . This 
suggests that LAEs have on average lower Ah than 
LBGs. 

For a lower Ah , Lya photons have a reduced chance of 
absorption by dust before escaping the galaxy because of 
a smaller number of resonant scattering. Thus, we infer 
that LAEs have strong Lya emission because of low Hi 
column densities. Although our data tell nothing about 
why LAEs have low column densities, possible reasons 
would in clude a high gas tempe rature caused by a shock 
heating (|Nakaiima et al . 2012b), which reduces the frac- 
tion of neutral hy drogen, and a clumpy gas geometry 
(|Bond et aLl[2O10f) . 



6. CONCLUSIONS 
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Fig. 12. — AuLya as a function of SFR (left panel), ff v (Ha) (middle), and £sfr (right). For all panels, the red and blue circles indicate, 
respe ctively, our LAEs and LBGs taken from the literature. The red tr iangles in the left pane l denote LAEs at z ~ 2.3 (Finkclstein ct al. 
[2011]), and the red squares in the middle panel indicate LAEs at z ~ 3 (McLinden ct al. 2011). 
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TABLE 7 

Summary of the physical quantities 



Object 


E{B - V), 


M« 


EW(Lyct) 




L(Hq) 


(T v (Ha) 


SFR(Ha) 


Jesc 


r 


'/ 






1O 9 M 


A 


km s _1 


10 41 erg s- 1 


km s — 1 


M Q yr" 1 




kpc 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




CDFS-3865 


n i«c;+0-009 
u.±oo_ 009 


n iq+0.21 
d - 18 -0.13 




190+H 


136.5+gf 


110.8 


190+13 


n 14+0.03 
u ' i4 -0.03 


1 1+ ' 1 
A - 0.1 




CDFS-6482 


n is^+ - 026 
o.i«a_ 01g 


5 30 +1 18 






34.51^ 


67.5 


48+ 10 


29+ ' 20 
u - zy -0.20 


°- u -0.8 




COSMOS-13636 


n 97 o+0.018 

u.z(j_ 07g 


1 gq+0.39 


73+5 




9-5±i;l 


30.7 


18+| 


"57+ - 13 
u -°'-0.13 


Q n+1-8 
9 - u -4.8 


1Q+ 007 


COSMOS-30679 


n c9s+ - 0026 
u.ozo_ 0026 


19.75±|;|3 


87± 7 7 






28.1 


4 5 +_l 


i 7 +0.03 
u - 1 '-0.03 


107.5±|?;« 


o so+°- 03 



Notes. (1) Dust extinction estimated from SED fitting; (2) Stellar mass estimated from SED fitting; (3) Rest-frame Lya EW; (4) Velocity offset 
of the Lyct line; (5) Ho luminosity; ( 6) Ve locity dispersion of Ha line; (7) S tar f ormation rate from Ha luminosity corrected for dust extinction; 
(8) Lyo escape fraction estimated in q4.5l (9) Size of the g alaxy derived in i|4. 61 under the assumption that stars dominate the total mass in the 
luminous part of the galaxy; (10) dumpiness parameter f £|5,3.2ll . 
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We have presented the results of Magcllan/MMIRS 
and Kcck/NIRSPEC spectroscopy for five LAEs at z ~ 
2.2 for which high- resolution FUV spectra from Magcl- 
lan/MagE are available. These objects are taken from 
the z ~ 2.2 LAE samples constructed by K. Na ka.jima 
et al. (in preparation: see iNakaiima et al.l[2012al for the 
selection). The redshift of z — 2.2 is unique since we can 
observe from the ground all of Lya, LIS absorption lines, 
and optical nebular emission lines including Ha. 

We have successfully detected Ha emission from all 
five objects, and [On] AA3726,3729, H/?, and/or [Om] 
AA4959, 5007 for some LAEs on the individual basis. In 
addition to that, we have detected LIS absorption lines 
in a stacked FUV spectrum. After removing an AGN- 
contaminated object, we have measured the velocity off- 
sets of the Lya line (A^Lya) and of LIS absorption lines 
(A?j a bs) from the systemic redshift determined by nebu- 
lar emission lines, to discuss the gas motions of LAEs. 
The major results of our study are summarized below. 

• We have obtained Av Lya = 190±?| km s" 1 (CDFS- 
3865), 75±|§ km s" 1 (CDFS-6482), 18 ± 16 km 
s- 1 (COSMOS-13636), and 170 ± 16 km s" 1 
(COSMOS-30679); three out of the four have a 
positive Ai>L ya beyond the 2a uncertainty. Com- 
bining with the result that all four have a posi- 
tive weighted skewness, we have conclude that the 
Lya emission of our objects is mostly originated 
from outflowing gas. This conclusion is supported 
by the finding that the LIS absorption lines in 
the stacked FUV spectrum are blue-shifted with 
Au abs = -179 ± 73 km s _1 . 

• For a sample of eight z ~ 2-3 LAEs without AGN 
from our study and the literature, we have obtained 
Az^Lya = 175 ± 35 km s _1 , which is significantly 
smaller than that of LBGs, Ai>L ya — 400 km s" 1 . 
If LAEs at z > 6 have similarly small Avh ya values, 
some reionization studies based on LAEs assuming 
A^Lya as large as that of LBGs may be overesti- 
mating the neutral fraction of the IGM. 

• We have found an anti-correlation between 
EW(Lya) and AvL ya in a compilation of LAE, 
LAB, and LBG samples, i.e., EW(Lya) decreases 
with Ai>Lya- On a reasonable assumption that 
Ai>L y a positively correlates with the outflow ve- 
locity within a factor of two, this anti-correlation 
indicates that high outflow velocities are not the 
physical origin of the strong Lya emission of LAEs. 

• We have found that LAEs have comparable to or 
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smaller outflow velocities than LBGs at a given 
SFR, a(Ha), and £sfr (when the systematic er- 
ror of < 2 in AuLya as a measure of u out is con- 
sidered). We have also found that the slope of the 
AuLya-SFR relation is different between LAEs and 
LBGs. Thus, the physical origin of LAEs' outflows 
appears to be different from that of LBGs'. It ap- 
pears that LAEs' outflows are at least consistent 
with the radiation pressure case. 

• To identify the physical origin of large-EW(Lya) 
galaxies, we have tested three hypotheses which 
may facilitate the escape of Lya photons: weak 
dust extinction, clumpy-cloud geometry, and low 
Ah- Since we have found no significant correla- 
tion between EW(Lya) and E(B — V), it is not 
clear if weak dust extinction leads to a high escape 
fraction of Lya photons. Although there may exist 
an anti-correlation between f^/ c a and the dumpi- 
ness parameter, more data and theoretical work are 
needed to draw a firm conclusion. 

• We have found that LAEs have Av-L ya ~ Av a b s |, 
in contrast with LBGs which have A^Lya ~ 2— 
3 x |Aw a b s |. When combined with the simulations 
of Lya radiative transfer i n a galaxy with an ex- 
pand ing shell of the ISM bv lVerhamme et"aT1l2006L 
I2008L this suggests that the typical A H of LAEs 
would be lower than that of LBGs, giving a smaller 
number of resonant scattering. Such low Ah may 
cause the observed strong Lya emission of LAEs. 
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